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Abstract 

The Green function of the fractional Laplacian of the differential 
order bigger than one and the Green function of its gradient perturba- 
tions are comparable for bounded smooth multidimensional open sets 
if the drift function is in an appropriate Kato class. 

1 Introduction 

Perturbations of the Laplace operator A by the first order or gradient op- 
erators b(x) ■ V were studied by Cranston and Zhao in [23]. They proved 
for Lipschitz domains that the Green function and the harmonic measure 
of A + b(x) ■ V are comparable with those of A under an appropriate Kato 
condition on the drift function b. Zhang then showed in [16] and (17] that 
the transition density of A + b ■ V has Gaussian bounds. The results were 
extended to more general second order elliptic operators by Liskevich and 
Zhang Q40J), and to drift measures satisfying the Kato condition by Kim 
and Song ([36]). 
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The fractional Laplacian A Q / 2 , < a < 2, is a primary example of a non- 
local generator of a Markovian semigroup. Perturbations of A a / 2 received 
much attention recently. In particular Schrodinger perturbations of A"/ 2 
were studied by Chen and Song ( [19] , [21] ) , Bogdan and Byczkowski ( [7] , [8] ) , 
Bogdan, Hansen and Jakubowski ([TT]) and Bogdan, Burdzy and Chen ([6]). 
Non-local Schrodinger-type perturbations were considered by Kim and Lee 
in [55], following earlier papers of Song ([H]. [IS]). Gradient perturbations 
of A 1 / 2 were studied by Caffarelli and Vasseur ([16]) and Kiselev, Nazarov, 
Volberg ([37]). Gradient perturbations of A Q//2 for a > 1 were considered 
by Bogdan and Jakubowski ([12]) and Jakubowski and Szczypkowski ([31]). 
with focus on sharp estimates of the corresponding transition densities on 
the whole of M d . In the present paper we estimate the Green function for 
smooth bounded subsets of M. d . 

Following [12] we let a G (1, 2). We will consider dimensions d G {2, 3, . . .}, 
a nonempty bounded open C 1 ' 1 set D C W 1 , its Green function Gd for A Q//2 , 
and the Green function Gd of the operator 

L = A Q / 2 + b{x) ■ V , 

where b is a function in Kato class JC^~ (for details see Section [2]). Our 
interest in L is motivated by the development of the classical theory of the 
Laplacian, non-symmetry of L (we have L* = A a / 2 — b(x) ■ V — divfr), the 
fact that the drift is quite a problematic addition to a jump type process, 
and by a handful of techniques which already exist for A"/ 2 . 

The following estimate, aforementioned in the Abstract, is an extension 
to A a / 2 of the results of Cranston and Zhao |23j . 

Theorem 1. Let d > 2, 1 < a < 2, b G K, a d ~ x , and let D C R d be bounded 
and C 1 ' 1 . There exists a constant C = C(a, b, D) such that for x,y G D, 

C-'Gnix, y) < G D (x, y) < CG D (x, y) . (1) 

Sharp explicit estimates of Gd, hence of Gd, exist, see §M§, and sharp 
explicit estimates of the corresponding Poisson kernel are given in f!72|) below. 
Theorem [1] is based on the perturbation formula for the Green operators, 

Gd — Gd + Gd frV Gd , 
where feV tp{x) = b(x) ■ Vip(x). Iterating yields formal perturbation series, 

oo 

G D = J2G D (bVG D ) n . 

n=0 
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The structure of the proof of Theorem [T] is now as follows. Section [2] provides 
details on the C 1 ' 1 condition and on transition densities, Green kernels and 
harmonic functions of the underlying Markov processes. In Section[3]we prove 
the perturbation formula and in Section H] we prove that the perturbation 
series indeed converge and yield ([T|) for small sets D with bounded distortion. 
In the proofs we use estimates for Go ([3B], [2D], [SO]) and for the gradient of 
Go (|14j). the boundary Harnack inequality for A Q / 2 ([13], [5]) an d the Kato 
condition (|30|) for the drift function b. As a result in Section |4] we obtain 
the Harnack and boundary Harnack inequalities for nonnegative harmonic 
functions of L in large open sets. These are then used in Section [5] along 
with the perturbation formula and a rough upper bound for Go given in 
Lemma to prove Theorem [T] for arbitrary bounded C 1 ' 1 open sets. A 
number of other auxiliary results are proved in the Appendix. 

Concerning the statement of Theorem [H we note that if the diameter of D 
is smaller than r and the distortion of D is smaller than A, then the constant 
C in ([1]) depends only on d, a, r, A and the suprema in the definition of the 
Kato class /Q (see below). 

We observe that an approach similar to ours was recently used for gradi- 
ent perturbations of elliptic operators on small sets in [28] (see also [16]). In a 
wider perspective, Theorem [1] is an analogue of the Conditional Gauge Theo- 
rem (CGT) in the theory of Schrodinger perturbations, see [TJ, [8], [19], [21], 
[22] and [25]. We should remark here that the distributions of the Markov 
processes generated by A a / 2 and L are not mutually absolutely continuous 
locally in time even for (nonzero) constant drift b, and any a G (0,2), see 
[121 Theorem 33.1]. Therefore techniques based on the Girsanov theorem 
([23] ) seem unavailable, and we need to proceed via analytic estimates of 
kernel functions. Apparently an adaptation of our arguments could be used 
to give a short analytic proof of CGT (compare [22], [Z], P23), i n fact a proof 
much simpler than that of Theorem [IJ Noteworthy, Green function estimates 
for Schrodinger perturbations hold conditionally under global assumptions 
of finiteness, e.g. gaugeability, existence of (finite) superharmonic functions 
bounded from below or smallness of the spectral radius. Lemma [7J, a con- 
sequence of the estimates of the transition densities in [12], overrides such 
assumptions here. Heuristically, adding drift b(X t )dt to a stochastic process 
will not increase its mass on M. d . This contrasts with a possibly exponential 
growth of the mass of Feynman-Kac semigroups generated by Schrodinger 
operators. The drift may, however, change the mass of the process killed off 
D by trying to push it away from the fatal D c . This is why ([T]) is nontrivial 
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phenomenologically. Also the symmetry of the semigroup and Green function 
are lost in the presence of the drift, causing certain technical problems. In 
this connection we note that Gd{u, %) may be considered the Green function 
of L*, and this operator has non-zero Schrodinger part, namely — divfr. Our 
results apply in particular to the Ornstein-Uhlenbeck operator A"/ 2 + kx ■ V 
(for dimensions d > 2 and 1 < a < 2). Here k is a constant. We refer to [21] 
and [32J for estimates of superharmonic functions of this important operator. 
We note that for k < the drift function b(x) = kx will generally increase 
the occupation time density (i.e. the Green function) for sets D containing 
the origin. 

The proof of (|T]) turned out to be quite difficult to handle, in terms 
of both the preliminaries and the auxiliary estimates of the Green function. 
Therefore we focused our attention on the more explicit C 1 ' open sets rather 
than Lipschitz open sets. We hope that our approach may now be adapted 
in the Lipschitz case. Here the sensitive elements are Lemma |9] and (|47|) . 

A few additional comments on possible extensions of the results are due. 
If d — 1 < a, then then the right hand side of f[2"9"j) below will no longer 
be integrable. This explains our restriction to d > 2. We however con- 
jecture that Theorem [TJ does extend to d = 1. This case is interesting 
even for the sake of the one-dimensional Ornstein-Uhlenbeck process. One 
may wonder if (JTJ) holds for a = 1, but we certainly know that ([1]) fails for 
a G (0, 1). Indeed, if < a < 1 then the expected exit times from balls, 
to wit, f B ^ xor ^GB(x ,r)(x,y)dy, are generally incomparable for A Q / 2 and the 

Ornstein-Uhlenbeck operator A Q//2 + kx ■ V (see [SI]), so the Green functions 
are not comparable either. Heuristically, a (first-order) gradient perturbation 
is infinitesimally small with respect to A Q / 2 only if a > 1. This explains the 
restriction 1 < a < 2 in [12], [M] and the present paper. We remark that 
the existence of ratios and Martin representation of nonnegative harmonic 
functions of L may likely be obtained with the results and toolbox presented 
in this paper and [T3]. We also note that a similar approach should apply to 
additive perturbations of A"/ 2 by non-local Levy-type operators (compare 
[24]). provided ( 129]) can be generalized. It also seems possible and interesting 
to study drift perturbations of more general semigroups subordinated to the 
Gaussian semigroup ([4"3]). 
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2 Preliminaries 

In what follows, lR d denotes the Euclidean space of dimension d > 2, dy 
stands for the Lebesgue measure on M. d , and we let 

1 < a < 2. 

Without further mention we will only consider Borelian sets, measures and 
functions in M d . By x-y we denote the Euclidean scalar product of x, y G M. d . 
We let B(x, r) = {y G M. d : \x — y\ < r}. For D C M d we denote 

Sd(x) = dist(x, D c ) , 

the distance to the complement of D. 

Definition 1. Nonempty open D C lR d is of class C 1,1 at scale r > if for 
every Q G <9.D there are balls B(x',r) C -D and B(x",r) C -D c tangent at Q. 

Thus, B(x',r) and B(x",r) are the inner an cmter balls tangent at Q, 
respectively. If -D is C 1 ' 1 at some unspecified scale (hence also at all smaller 
scales), then we simply say D is C 1,1 . The localization radius, 

r o = r o(D) = sup{r : D is C ' at scale r}, 

refers to the local geometry of D, while the diameter, 

diam(D) = sup{|x — y\ : x,y & D} , 

refers to the global geometry of -D. The ratio diam(_D)/r (-D) > 2 will be 
called the distortion of D. We can localize each C ' open set as follows. 

Lemma 1. There exists k > snc/i t/jat if D is C 1,1 at sca/e r ana 1 Q G 
then there is a C 1 ' 1 domain F C D with ro(F) > nr, diam(F) < 2r and 

D n B(Q, r/4) = F n fi(Q, r/4) . (2) 

We will write F = F(z,r), and we note that the distortion of F is at most 
2/k, an absolute constant. The proof of Lemma [T] is given in the Appendix. 

In what follows D will be a nonempty bounded C 1 ' 1 open set in M d . 

We note that such D may be disconnected but then it may only have a finite 
number of connected components, at a positive distance from each other. 
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We will now give a brief review of the potential theory of the fractional 
Laplacian, and of the fractional Laplacian perturbed by gradient operators. 
The former case is well known ( [39] , [2] , [9] , [7J , [13] ) . The latter case is similar 
but we feel it calls for more details, and they are given in the Appendix. 

Let A da = r((d- 7 )/2)/(2T7r d / 2 |r( 7 /2)|) and 

u(y) = A dt . a \y\~ d ' a , yeR d . 
The coefficient Ad,- a is so chosen that 

[ [l-co8(Z'y)]u(y)dy=\Z\ a , £eR d . (3) 

For (smooth compactly supported) G C%°(M. d ), the fractional Laplacian is 

A a/2 <p(x) = lim / [<j){x + y)- <f>{x)\ v{y)dy , x e R d (4) 

£ +° J\y\>e 

(see [9] for a broader setup). If x ^ supp (p then 

A a ' 2 ^x) = / <j>{y)v{y-x)dy. (5) 

If r > and (f) r (x) = (j)(rx) then 

A a/2 r (x) = r a A a/2 <p(rx) , xeR d . (6) 

In this respect, A a ^ 2 behaves like differentiation of order a. We let pt be the 
smooth real- valued function on M. d with Fourier transform 

/ Pt (x)e ix < dx = e~ ma , t > , C e K d . (7) 

According to ([3]) and the Levy-Khinchine formula, {p t } is a probabilistic 
convolution semigroup with Levy measure u(y)dy, see [12], [15] or [9]. Let 

p(t,x,y) = ptiy - x) . 

Using ([7]) one proves that p is the heat kernel of the fractional Laplacian: 

oo 

//p( U - W )[^(M+A?'V(M] **.= -#(.,.), (8) 
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where s G R, x G R d and G C™(R x R d ). 

We consider the time-homogeneous transition probability 

(t, x, A) i-)- / p(t, a, , i>0,i6i d ,/lcR d . 

By Kolmogorov's and Dinkin-Kinney's theorems the transition probability 
defines in the usual way Markov probability measures {P*, x G R d } on the 
space Q of the right-continuous and left-limited functions u : [0, oo) — > R d . 
We let W be the corresponding integrations. We will denote by X = {X t } t > 
the canonical process on Q, X t (cu) = u(t). In particular, according to (J7J), 

E°e lXH = e~ m \ i e R d , t > . (9) 

In fact, (X,F°) is a Levy process in R d with zero Gaussian part and drift, 
and with u(y)dy as the Levy measure (12]. It follows from ([7]) that 

Pt ( x ) = r d/a Pl {r 1/a x) , t>o, xeR d . (10) 

It is well-known that pi(x) ~ 1 A |a;| _d_Q , hence 

^)^ rd/ttA ^P' i>0,i6l d . (11) 

^* 

Symbol means that either ratio of the sides is bounded by C G (0, oo), and 
C does not depend on the variables shown, here t and x. We will write mere 
~ if C is unimportant or understood. Constants will usually be denoted 
with generic C (in statements) or c (in proofs), and we will occasionally 
enumerate them for convenience of referencing. As usual, a A b = min(a, b) 
and a V b = max(a, b). In what follows we will often use the identity 

ab= (oA6)(oV6). (12) 

In view of (TlOl) and the fact that each p 4 is a radial function, X is called 
the isotropic ct-stable Levy process (see [12], [15] for a discussion of general 
stable Levy processes). We introduce the Riesz potential kernel (for d > a), 

POO 

A d , a \x\ a ~ d = / Pt(x)dt, xeR d . (13) 
Jo 

This is infinite if x = 0, see (fTTj) . 
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To study A"/ 2 with Dirichlet conditions we will consider the time of the 
first exit of the (canonical) process from D, 

r D = inf{t > : X t i D} . 

We let u x D {B) = P X {X TD G B), the a-harmonic measure of D (0, g], gH]). 
The joint distribution of (td,X Td ) defines the transition density of the pro- 
cess killed when leaving D ([27], [3], [22]): 

p D (t,x,y) =p(t,x,y)-E x [r D < t; p(t - r D: X TD ,y)\, t > 0, x, y G R d . 

By Blumenthal's 0-1 law, radial symmetry of pt and C 1 ' 1 geometry of the 
boundary of dD, we have P x (t£> = 0) = 1 for every x G -D c . In particular, 
Pn(t, x, y) = if x G D c or y G -D c . By the strong Markov property, 

E x [t < m; f(X t )} = I f(y)p D (t,x,y)dy, t>0,xGM d , 

for functions / > 0. The Chapman-Kolmogorov equations hold for p^, 

Pd(s, x, z)p D (t, z, y)dz = p D (s + t, x, y) , s, t > 0, x, y G M d . 



Also, pd is jointly continuous when t ^ 0, and we have 

< x, y) = pz)(t, y, x) < p(t, x, y) . (14) 

In particular, 

/ PD(t,x,y)dy <1. (15) 

</R d 

For s G R, x G M d , and G C£°(M x D), we have (compare (JSJ)) 

oo 

//p D («-.,.,,)M**) + A?»««,,)] ^=-^,,), (16) 

which justifies calling the heat kernel of the (Dirichlet) fractional Lapla- 
cian on D. We define 



oo 







G D (x,y)= I p D (t,x,y)dt, x,yeR d . (17) 
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It follows that Gu(x,y) is symmetric and lower semi-continuous, and 



G D (x,y)+ / A d>a \y - z\ a d u x D (dz) = A d>a \x - y 



I a— d 



(18) 



D c 



The Green operator of A Q / 2 for D is 



TD 



f(X t )dt 



G D f(x 
and we have 

G D (A a/2 <p){x) = -cj)(x) , xeR d 



GD(x,y)f(y)dy, xe 



e C?{D) 



(19) 



A result of Ikeda and Watanabe [29] asserts that for x G D the P x -distribution 
of (td, X Td -, X Td ) restricted to X TD _ ^ X TD is given by the density function 



(s, u, z) !->■ x, u)u(z — u) . 



(20) 



The C 1 ' 1 geometry of D implies that P X (X TD ^ ^ X TD ) = 1 for z e D (0). 
By (fT7|) . (1201 and Tonelli's theorem the P x -distribution of X TD has a density 
function, called the Poisson kernel and defined as 

P D (x,z)= [ G D (x,y)is(z-y)dy. (21) 

J D 

The Green function and Poisson kernel of the ball are known explicitly: 

G B (x ,r)(x,v) = B d)a \x - v\ a ~ C 



w s a/2-l 



s + iy/2 as ' 



PB(x ,r){x,y) = C dja 



r — \x — xo\ 
\y — Xq\ 2 — r 2 



a/2 



\x-y\ 



(22) 
(23) 



where B dyCl = T(d/2)/(2 a 7r d / 2 [T(a/2)] 2 ), C d , a = T(d /2)7r~ 1 ~ d / 2 sm(vra/2), 

w = (r 2 — \x — xo| 2 )(r 2 — \v — xo\ 2 )/\x — v\ 2 , 

\x — Xq\ < r, \v — Xq\ < r, and \y — Xq\ > r; see [1], jH] or [39] . 

The next estimate was proved by Kulczycki [38] and Chen and Song [20] . 



Go(x,y) 



c 



\x-y 



a—d 



5 D {xY' 2 5 D {yrl 2 



A 1 



f - y\ 



I a— <i 



f — y\ 

[5 D (x) V \x-y\ V5 D (y)] a ' 
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x,y E D . 



(24) 



(25) 



The reader may check equivalence of (1241) and (I25p by first considering the 

3 

case 5d(x) ~ So(y)- We like to remark that f[2"o"j) may be also regarded a 
direct consequence of the approximate factorization of the Green function 
of Lipschitz open sets, see [30j Theorem 21]. It is well known that C = 
C(d, a, A) in ([21]) and fl25]), if diam(D) / r (D) < A, i.e. C may be so selected 
to depend only on d, a and (an upper bound for) the distortion of D. This 
follows from the proofs of [38] and [2U] and is explicitly stated in [30J, see 
also [35] . 

We will consider a nonnegative function m on R d , and an open set U cM. d . 
u is called a-harmonic on U if for each open bounded V C V C U, 

u(x) = E x u(X Tv ), xeV. 

We say that u is regular a-harmonic on U if also 

u(x) = Wu(X Tu ) , x eU . 

Here we assume absolute integrability of the expectations, and K x u(X T[/ ) 
is understood as E x [t[/ < oo; u(X Tu )]. For instance x \-> Go(x,y) is en- 
harmonic in D \ {y}. In fact, by the strong Markov property, Gu(x,y) = 
Gy(x,y) + WGd{X tv , y) for every open V C U, and Gy{x,y) = if 
dist(y,F)>0(see, e.g., [13]). 

The following two results can be found in [5], see also [T3] . 

Lemma 2 (Harnack inequality) . Let x, y G M. d , s > and k G N satisfy 
\x — y\ — 2 fc s. Let function u be nonnegative in M. d and a-harmonic in 
B(x,s) UB(y,s). There is C = C(d,a) such that 

C -^2' k ^ u (x) < u{y) < C2 k{d+a) u(x). (26) 

Lemma 3. Let Z G dD andr G (0, ro], < p < 1. Assume that functions u, 
v are nonnegative in M. d and regular a-harmonic and non-zero in DC\B(xq, r) . 
If u and v vanish on D c D B(xo,r) then 

c ^u{x) M 

v{x) v[y) v{x) 

for x,y G D PI B(xo,pr). Here C = C(d, a,p). 
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We like to remark that the boundary Harnack inequality (Lemma [3]) in 
fact holds for general open sets and is equivalent to an approximate factor- 
ization of the Poisson kernel of general open sets, see [13]. We encourage 
the reader to factorize Pb(o,i)(x, y) when x, y are not too close to each other. 
In passing we also note that an approximate factorization oi po(t,x,y) for 
Lipschitz domains is given in [ID] . Concluding this part of our preliminary 
discussion we refer the reader to [5], [13] for more details and references. 

We note that a-harmonic functions are smooth where a-harmonic; use 
( 123]) or see [7]. The following gradient estimate is given in [T4"| Lemma 3.2]. 

Lemma 4. Let U be an arbitrary open set in lR d . For every nonnegative 
function u on M. d which is a -harmonic in U we have 

\Vu(x)\ < dpj\ , xeU. (28) 

Since Gu(-, y) is a-harmonic in U \ {y}, for every y G U we obtain 

\V x Gu(x,y)\<d G "[ X ; y) x,yeU,x^y. (29) 

djj{x) A \x -y\ 

We note in passing that a reverse inequality holds locally at the boundary of 
Lipschitz domains, with constant depending on the Lipschitz character of D 
( [T4"l Lemma 4.5]). In this sense (|28|) and (129]) are sharp. Also, V x Gu(x,y) 
is jointly continuous for x ^ y G U, see [TH (10)]. 

Recall that 1 < a < 2. We say that vector field b : R d — > M. d belongs to 
the Kato class /Q -1 if 

limsup / \b{z)\ \x-z\ a - l ~ d dz = 0. (30) 

£ ^°xm d J\x-z\<e 

For instance, if b is bounded or if \b(z)\ < and < e < a — 1, then 

b G K%~ . Without much mention elements of /Q -1 will either be vector 
fields R d — > M. d or real-valued test functions R d — > R, i.e. is more a 

condition than a class. Since \x — is locally bounded from below, 

\b(z)\dz is a locally finite measure, and (1301) is a local uniform integrability 
condition. If b G /C^ -1 and / is bounded, then fb G in particular, fb 

is locally integrable. We note that /C^ _1 C tC°[, where /Q is defined by 

limsup/ \b(z)\ \x - z\ a ~ d dz = . (31) 

£ ^°x£R d J\x-z\<£ 
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Following [12] and [M] we recursively define, for t > and 1,1/6 M. d , 



Po(t,x,y) =p(t,x,y) , 
p n (t,x,y)= / p n -i(* - 2)6(2) • V z p(s,z,y)dzds, n>l, 

Jo JR d 



and we let 

oo 

P=Y,Pn- (32) 
n=0 

The series converges absolutely, p is a continuous probability transition den- 
sity function, and 

c^p(t, x, y) < p(t, x, y) < crp(t, x, y) , x, y G R d , < t < T , (33) 

where c-r — > 1 if T — > 0, see [121 Theorem 2]. From a general perspective 
the approach of [12], [31] consist of using the semigroup as test functions, 
setting the assumptions on the perturbation so that p\ is dominated by p in 
short time, and recursively estimating multiple integrals defining p n , so that 
the comparability with p is preserved. Auxiliary estimates of V 2 p(s, z, y) are 
obtained in [T2"l [31] by subordination to the Gaussian kernel, but the scope of 
the method is wider. For instance applications to Schrodinger perturbations 
of general transition densities are given in [11], [33] . 

We let P, E be the Markov distributions and expectations defined by 
transition density p on the canonical path space. We define the heat kernel 
of L on D by the usual G. Hunt's formula, 

p D (t, x, y) = pit, x, y) - W [r D < t; p(t - r D , X TD , y)] . (34) 

We denote by Gd(x, y) and Gd the Green function and operator of L on D, 

POO 

G D (x,y)= p D (t,x,y)dt, (35) 
Jo 

G D 4>(x)= G D (x,y)(f)(y)dy. 

By Blumenthal's 0-1 law, pi)(t,x,y) = and Go{x,y) = if x G D c or 
y G D c , see ( )33|) . The next lemmas rely on the definition of p and generalize 
results stated above for A Q / 2 . The proofs of Lemma [5], [6] and [8] are moved to 
the Appendix. 
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Lemma 5. For s > 0, x G D and <f> G C£°((0, oo) x D) we have 
J J Pd(u-s,x,z) (d u + A* /2 + b(z) ■ V z ) (f) (u, z) dz du = -<p(s,x) . (36) 

By fl33|) we have 

r P(t,x,y) P(t,x,y) , . 

hm = hm = v [y — x . 

i^o t t-»o t 

Thus the intensity of jumps of the canonical process X under is the same 
as under F x . Accordingly, we obtain the following description. 

Lemma 6. The F x -distribution of (td, X td ) on (0, oo) x (D) c has density 

/ p D (u,x,y)v(z -y)dy , u>0,5 D {z)>0. (37) 
Jd 

We define the Poisson kernel of D for L, 

P D (x,y) = / G D (x,z)v(y - z) dz , x G D , y G D c . (38) 

JD 

By ((35D, (EHD and (EZD we have 

¥ X (X TD EA)= [ [ G D {x,z)v{y- z)dzdy = [ P D (x,y)dy, (39) 
JaJd J a 



if A C (D) c . For the case of A C <9-D, we refer the reader to Lemma [T] 

The following rough estimate of Go results from the estimates of p and 
the fact that X jumps out of D at least with intensity JLi >diam rD) u (y)dy > 0. 

Lemma 7. Go(x,y) is continuous for x ^ y, Gd{x,x) = oo for x G D, and 

G D (x,y)<C Q \x-y\ a - d , x,yeR d , 

where Co = Co(d, a, diam(D)). 

Proof. We claim that there are constants c and C such that 

PD(t,x,y) < Ce~ c \ t>l, x,yeR d . (40) 
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Indeed, let kd(v) = J D c v{z — y)dz, so that kd(u) > c > for y G D. Let 
x e D, t > 0, and F(t) = f x \r D > t) = J D p D (t, x,y)dy. By Lemma El 



-F'(t) = [ ; 

Jd 



p D (t, x, y)K D (y)dy > cF(t) , 

Jd 

hence ¥ x {jd > t) < e~ ct . By the semigroup property and (13"B"|) . for t > 1, 

pD(t,x,y)< / p D (t-l,x,z)p(l,z,y)dz 
Jd 

< ci / p D {t-l,x,z)p(l,z,y)dz 



Jd 

< c lP {l, 0, 0) P x (r D > t - 1) < Ce" ci 



By (|35|). ([33]) and ([40]) we obtain 

/l /*OG 
c 1 p(t,x,y)dt + / Ce~ ct dt 

< A d , a \x - y\ a ~ d + C/c < (A*,a + C diam(D) d ~ a /c) \x - y\ a ~ d . 

By (!35|) and dominated convergence theorem Gr>(x, y) is continuous if x 7^ y, 
see (El and ([II]). □ 



The next lemma results from integrating (]36l) against time. 
Lemma 8. For a// ip e C^°(D) and x £ D we have 

G D {x, z) (A a / 2 v (z) + 6(z) ■ Vp(z)) cfe = -cp{x) . (41) 



The definition of L-harmonicity is analogous to that of a-harmonicity. 

Definition 2. u is L-harmonic on U if for each open bounded V dV dU , 

u(x) = Wu{X TV ), x 6 V. 

We say that u is regular L-harmonic on U if also 

u(x) =E x u(X TU ) , xeU. 

Here E, x u(X TU ) = E, x [tjj < 00; u(X TU )] and we always assume absolute 
integrability. In particular, x 1— > Gr>(x,y) is L-harmonic in D \ {y}, in fact 
Gd(x,v) = Gu(x,y) + WG v (X(Tu),y) for every open U C D. We should 
note that in general Gn(x, y) 7^ Go(y, x) (non-symmetry), and y \-> Gc(x, y) 
is not L-harmonic. This accounts in part for the difficulties in estimating Gr>. 
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3 Perturbation formula 

As before, D is a bounded C ' open set in R d , b G and 1 < a < 2 < d. 
Let 

G = Gd, G = Gd and 5 = 5d- 

In view of ( 1291) the next lemma yields uniform integrability of b(z) ■ V z G(y, z). 
In particular, the singularity 5(z) a l 2 ~ l of V Z G at dD integrates against \b\. 

Lemma 9. G(y, z)/[5(z) A|y— z\] is uniformly in y integrable against \b(z)\dz. 

Proof. In view of (125 p it is enough to prove the uniform integrability of 



H(y, z) = \y-z 



a _ d 5{y)^5{z)^ 1 



[5(z)V\y-z\v5(y)^5(z)A\y-z\ ' 
Let A R (y) = {z G D: H(y, z) > R} for R > 0. We will verify that 



lim sup / H(y, z)\b(z)\ dz = . 
?->oc eD J A , > 



For r > we denote 

K r = sup / IK^Hx-yr" 1 -^. 

xGR d JB(x,r) 

By ([30]) we have that < oo and K r 1 as r 1 0. For all x G M ' and r > 0, 
|6(z)| dz < r d+1 ~ Q / \x- z\ a ^~ d \b(z)\dz < K r r d+1 ~ a . 

B(x,r) JB(x,r) 

For r > we let D(r) = {z G D : 5(z) > r}. If m > 0, y G D and 
z G D(5(y)/m), that is <5(z) > 5(y)/m, then we have 

Indeed, by $12]), 

Q „ d %)«/ 2 ^) o/2 5(z)V\y-z\ 



H(y,z) = \y-z 



[6(z)V\y-z\V6(y)] a 5{z)\y - z\ 

< \y - z| a - 1 " ci 5(z) Q/2 - 1 <5(?/) Q/2 [5(z) V\y-z\V %)f~ Q 

< \y - z\ a - l - d 5(z) a/2 - l 5(y) l - a/2 < m l - a/2 \y - z^' 1 " 1 . 
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If R — > oo, then uniformly in y we have 
H(y,z)\b(z)\dz 

D(S(y)/m)nA R (y) 

<m x - a l 2 [ \y - z^^lbiz)] dz ^ . (42) 

J{zeD: \y-z\ a - 1 - d >Rm a / 2 - 1 } 

For y G D, k, n > and m > 2 we consider 
W™ k {y) = {zeD: -Mj < 5(z) < M, < | y _ ,| < (fc + i)^)}. 

may be covered by ci(A; + l) d ~ 2 m d ~ 1 2 n ( <i ~ 1 ) balls of radii thus 
/ |&(*)l<**< c 1 (ife+l) < '- 2 m d - 1 2 n ^^ sup / |6(z)| tZ« 

< c lK6(y)/m2n {k + l^nrt^ 

= c x K &{y)lm2n (k + l)^ 2 m^ 2 2"^ 2 )5(y) d+1 - Q . 

For 2 G W^ fc (?/) we have > 25(z), hence \y — z\ > 5{y)/2 and \y — z\ > 
S(z). We obtain 

/°° 00 /• "\a/2 

%«l<b<EE J \y-zMzY-^ )\ dz 

A R (y)\D(5(y)/m) n ~° k ~°W™ k (y) 

5{y) a ' 2 



^EE/ ; „,,,„X / ,,, / ^ Ml* 

n=0 fc=0 



- h (y) ((A; + l)<%)/2) a [%)/(m2«+i)]W2 

oo oo 

< c 2 K s{y)/m J> + l) d ^ Q-2 2 n(a - 2) 2 d (A; + i)- d 2-(W2) m i- a /2 

n=0 fc=0 

< ctfn af2 ~ x K S (y)i m ■ 

Let 5 > 0. We chose m and R so large that C3m a ' 2_1 i^ diani m)/ m < e/2 and 



sup / z)|6(z)| < e/2 . 

yeDjD(d(y)/m)nA R (y) 

This completes the proof. □ 
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We consider the operator oV : 

d 



bV(j))(x) = b(x) ■ V0(x) = J2 b > 



i=i 



dxi 



We will study the perturbation series Ylt^LoiGbV ) n G °f integral operators 
on /C^ _1 . Namely we will apply GbV G to real-valued / G /C^ _1 : 

GbVGf(x) = f G(x,z)b(z)-V z f G(z,y)f(y)dydz. (43) 
Jd Jd 

We will need to interchange the integration and differentiation in (|4"3"|) . 

Lemma 10. Let 1 < a < 2. If f E /C^ _1 or at least G(z,y)/\y — z\ is locally 
in z G D uniformly integrable against \f(y)\dy, then 

V z / G(z,y)f(y)dy = [ V z G(z,y)f(y)dy, zeD. (44) 
Jd Jd 

Proof. The result is proved for / G JC^ in [T4"l Lemma 5.2]. For the more 
general / we note that G(z, y) and V z G(z, y) are continuous on D x D 
except at z — y, see the remark following ([29]) . They are also uniformly 
integrable against \f(y)\dy for z in compact subsets of D. In consequence 
g(z) = Gf(z) and k(z) = f D 'V z G(z,y)f(y)dy are continuous on D. We 
consider f n = f A n V (— n) G fi'n = We have g n —> g and 

V(? n — » It follows that Vg = k. □ 

For x, ?/ G -D we let 

Jd G(x,y)(8(z)A\y-z\) 
u \ f \u m G(x, z)G(z, y)(5(x) A |x — yj) 

K{X,y) = / 6Uj 77; r777 \ i ix/r/ \ i rrCZZ. (46) 

v ' y; y D ' v G(x,y)(5(z) A\y - z\)(5(x) A\x - z\) 

In what follows k and k will serve as majorants for the perturbation series. 

Lemma 11. Let X,r < oo. There is C\ = Ci(d, a, b, A, r) swc/i £/ia£ i/ .D is 
C 1,1 , diam(D) /r (D) < A and diam(D) < r, then n(x,y) < C\, k{x,y) < 
2C\ for x,y G D, and Ci(d, a, b, A, r) — > as r — )• 0. 
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Proof. Denote g(x) = 5{x) a ^ 2 . Let x,z E D. By (I25p we have 

^l G(x z ) ~ 3\z) _ d 

g(x) ' (S(x) V \x-z\ V 



< r , . V r |x - z\ < —— — —- -\x-z\ a - d (47) 

~ \6(z) V \x - z\J - 8(z) V \x - z\ 1 

= (<f(z) A (x-^Dlx-^r" 1 ^. (48) 
Let Q(x,y) = G(x,y)/[g(x)g(y)}. The so-called 3G Theorem holds for Q: 

G(x,z) A j/) < cQ(x,y) , x,y,zeD, 
where c depends only on rf, a and the distortion of .D ([2S])- We obtain 

G(:n,*)G(z,y) - g2(z) ^f(^)< c fgM G(v) v 9M G (z,y) 



G{x,y) G(x,y) \g(x) g(y) 

t . 5(z) A \x - z\ 5{z) A\y — z\ 



\x-z\ d+1 - a \y-z\ d+l - a 
c(6(z) A \x - z\ A \y - z\) ( j- _ V |y _ ^ +1 „ a ) . (49) 



We have 

G(x, *)G(z, y) < c (|x _ zl"- 1 -" + | y - 2 |°-i-<*) 



G(x,y)(5(z) A \y — z\ A \x — z\) 

so we actually have uniform integrability against \b(z)\dz. The statement 
about k follows form ( |30|) . 

To estimate k we consider two cases. If 2\x — z\ > 8(x) A \x — y\, then 



S(z) A 


x — z 


A|y- 


- ^ 


1 


{5{z)A\y- 


-z\)(5(x)A 


x - 


- z\) - 8(x) A 


x — z\ 



S(x) A \x — y\ 

If 2\x — z\ < 5(x) A \x — y\, then 5(z) > 8(x)/2, \y — z\ > \x — y\/2, and so 

2 



6(z) A 


x — z 


A|l/- 


- z 


^ S(z) A \: 


TC — Z 


A\y-z 




(S(z) A \y - 


- z\)(5{x) A 


x - 


-z\)~ (5(x)/2A 


\x — 


vW\x- 


z\ 



< 



By fj4*9|) and fl30|) we obtain k(x,y) < 2C\. In fact we observe the uniform 
integrability against \b(z)\dz. The above estimates of the factors in (1451) and 
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(I4"6"j) depend on D only through d and diam(D)/r (^) ([3D])- Therefore the 
integrals in (I45p and (1461) are arbitrarily small if diam(D) is small enough, 
and the distortion of D is bounded by a constant. This follows from fl30|) . If 
diam(.D) is not small but finite then we only have the boundedness of k and 
k, which also follows from (I30p . □ 



For x 7^ y we let 

G 1 (x,y) = [ G{x,z)b{z)-V z G{z,y)dz. (50) 
Jd 

By (1291) . (145 jl and Lemma [TT] the integral is absolutely convergent, 

IG^yJI^dG^y) / 3^P^r^ dz < dC&fay) . (51) 

G(x,y)(d(z) A |y - z|) 



For / e K%~ 1 C /CJ we have 

< d / G(x,j/)|/(j/)|dj/<oo, 

hence by Lemma [TUl (T29|) and Fubini's theorem, 

GbVGf(x) = [ G(x,z) f b{z)-V z G{z,y)f{y)dydz 
Jd Jd 

= [ G 1 (x,y)f(y)dy. 
Jd 

We like to note that the linear map / t— > bVGf preserves /C^ _1 because VGf 
is a bounded function, see Lemma M and the remarks following (1HU|) . 
We will now prove the pointwise perturbation formula. 

Lemma 12. Let i,i/6 IR d ; x ^ y. We have 

G(x, y) = G(x, y)+ f G(x, z)b(z) ■ V z G{z, y)dz . (52) 
Jd 
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Proof. Let x G D. For (p G /C? we consider 



G(x,y) - G(x,y) - / G(x, z)b(z) ■ V z G(z, y)dz 



(j>{y)dy. 



By Lemma [9] and Lemma the iterated integral converges absolutely. If 
<p G C™{D) and <p = A a / 2 <p, then using (l4"Ij) . (115)1 and Lemma [TO] we obtain 

A(0) = - G(bV<p)(x) + p(a;) - G(6V(-<p))(x) = . 

By [3 Theorem 3.12], A(0) = L d <f)(y)X(y)dy, where A is a-harmonic on D. 
Since our A is bounded near 3D and vanishes on D c , we have that A = 0, see 
Lemma 17]. By uniform integrability and the remark following (1291) and 
by Lemma [7] we see that both sides of (1521 are continuous in y G D \ {x}, 
hence we have pointwise equality in (152|) . □ 

In addition to G\ we inductively define 
G n (x,y) = J G n -x(x,z)b(z) -V z G(z,y)dz , x^yeD, n = 2,3, .... 
We also let G (x, y) = G(x, y). By (J2"9"j) . (15T1) . (T4"5l) . LemmaHHand induction, 

^(C^)"" 1 / |^)|G(x,^)|V,G( 2 ,y)|^<(C lC /)"G(x,y), (53) 

./£> 

where n — 0, 1, 2, . . . and x ^ y. By (1551) and induction we prove that 



G n (x,z) / b(z) ■V z G(z,y)f(y)dydz = / G n+1 (x,y)f(y)dy, 

D J D J D 

hence G n is the integral kernel of G( bV G) n , 

G(bWG) n f(x)= [ G n (x,y)f(y)dy, f G /C^ 1 , x & D . 
Jd 

We can also handle the gradient of G n . Namely, for x,y G D, x ^ y, and 
n — 1, 2, . . ., we have: 



iV.Gn-iC*, < (2C 1 )"- 1 d w -. fjfj' y) , , (54) 

o(x) A |x — y\ 
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G n (x,y) = J G(x,z)b(z) ■V z G n -i(z,y)dz , (55) 

and 

V x G n (x, y) = J V x G n -i{x, z)b(z) ■ V z G{z, y)dz . (56) 

The inequality and the equalities are proved consecutively by induction. In 
the process we use Lemma [10], (|46|) . estimates following ([49]) in the proof of 
Lemma [TT1 and Fubini's theorem. 

Iterating (|52|) . by ([55]) we obtain for n = 0, 1, . . ., and x ^ y, 



G(x,y) = G(x, y) + J G(x, z)b(z) ■ V z G(z, y)dz 

n „ 

G k (x, y) + J G(x, z)b(z) ■ V z G n (z, y)dz . (57) 



fc=0 

The details are left to the reader. 



4 Local results 

We will prove and use the comparability of G$ and G$ for small smooth 
sets, S, of class C 1 ' 1 (this part of our development is similar to [28]). The 
following is a variant of Khasminski's lemma ([22]). 

Lemma 13. Let d > 2, 1 < a < 2, b e /Q -1 and A > 0. There is 
e = e(d,a,b, A) > such that if diam(S) /ro(S) < A and diam(S') < e, then 

^G s ( Xl y) < G s (x, y) < ^G 5 (x, y), Xj yeR d . (58) 

Proof. Let D = S. By Lemma dU and ( 153]) there is e = e(d, a, b, A) > and 

\G n (x,y)\ <4- n G(x,y), x ^ y , n = 0, 1, . . . , 

provided diam(D) /ro(D) < A and vq(D) < e. For x ^ y we have G(x,y) = 
Y^ =Q G n {x,y). Indeed, the remainder in ( 157|) is bounded by 




|6(z)|(2Ci) w -V r . r \ 4 z -)• , asn->oo. 
A |y - z| 
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Here 2C\d < 1/2 and the integral is finite because of Lemma and (13 ip . 
Thus, 

oo oo i 

G(x,y) = ^G n (x,y) < j^Gfoy) = g^y), 

n=0 n=0 

and 

oo 

G(x, y) > G(x, y)-J2 4 ~" G fo V) = a ' 

n=l d 

□ 

We like to note that the comparison constants in the above proof will 
improve to 1 if diam(S') — > and the distortion of S is bounded. By fl39|) . 

F X (X TS e A) « P X (X TS ei), ieS, Ac (S) c . (59) 

We are in a position to prove that the boundary of our general C 1 ' 1 open 
set D is not hit at the first exit (recall that 1 < a < 2). 

Lemma 14. For every x E D we have that F X (X TD £ <9-D) = 0. 

Proof. Let w(x) = P x (X rr , £ <9-D), x £ M d . We claim that there exists 
c = c(d, a, D, b) > such that u(x) < 1 — c for x £ -D. Indeed, we consider 
small £ > 0, x £ -D, r = e dist(x, _D C ), the ball B = B(x,r/2) C D, and a 
ball 5' C (D) c with radius and distance to B comparable with r. By ( 159]) 
and (1231) . 

P*(X TD £ <9£>) > P a: (X TB(a;)r . /2 ) £ B') « P*(X TB(a , r/2) £ S') > c. 

Furthermore, let _D n = {y £ D : dist(y,D c ) > 1/n}, n = 1,2,.... We 
consider n such that B(x, r/2) C D n . We have P x (X TDn £ D) < 1-¥ X (X TB £ 
B') < 1 — c, as before. Let C = sup{u(y) : y £ D}. We have u(x) = 
W{u{X TDn )- X TDn ED}< C(l - c), hence C < C(l - c) and so C = 0. □ 

In the context of Lemma [13], the P x distribution of X TS is absolutely 
continuous with respect to the Lebesgue measure, and has density function 

P s (x iy )^Ps(x,y), yeS c , (60) 

provided x £ S. This follows from (139]) and Lemma [TU For clarity, 

F X {X TS £ A) « P x (X rs £ A) , x £ 5 , A C S c . (61) 
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Lemma 15 (Harnack inequality for L). Let x,y G M. d , < s < 1 and 
fceN satisfy \x — y\ < 2 fc s. Lei tt fre nonnegative in M. d and L-harmonic in 
B(x, s) U B(y, s). There is C = C(d, a, b) such that 

c -x 2 -k( d+a )~^ < < C2 k{d+a) u(x) . (62) 

Proof. We may assume that s < 1 A e/2, with e of Lemma [131 Let tt(^) = 
for z G B(y,2s/3) c and u(z) = J B{y 2s/3) c u(v)P B ( y ,2s/3)(z, v) dv for z G 
B(y,2s/3), so that w is nonnegative in IR d and (regular) a-harmonic in 
B(y,2s/3). Let z G B(y,s/2). By 



u(*) = E 2 M(X(r B (j /i2s /3))) = / u(v)P B (y j2 s/3)(z,v) dv « li(z) . 

By Lemma [2] we get ~ with a constant depending only on d, a and 
b. To compare u(x) and ?/(?/) we will assume that \x — y\ > 3s/2, because 
otherwise we may take smaller s. For z G B(y,s/2) we have |x — z| < 
\ x — y\ + \y — A < 2 fc s + s/2 < 2 fc+1 s, and we get 

/ (s/2) 2 \ a/2 1 

PB {X ,S/2){X,Z) = C ^[ lx _ zl 2_ {s/2) 2) ^37p 

> C d , a 2- a s a \x- z\- d ~ a 

> C dya 2~ {d+2a) 2^ d+a) s~ d . 

Since Pb(x,s/2) ~ Pb(x,s/2), by the first part of the proof we obtain 

u(z) = / Pb(x,s/2)(x, z)u(z)dz > / P B{XtS/2 ){x,z)u{z)dz 

JB(x,s/2) c JB(y,s/2) 

« / P B{XiS/2 ){x,z)u{z)dz 

JB(y,s/2) 

> \B(y, s/2)\C d ^ d+2 ^2- k ^ d+a h- d u{y) = c2-^ d+ ^u{y) . 
By symmetry, £t(x) ~ □ 

We obtain a boundary Harnack principle for L and general C 1, sets Z). 

Lemma 16 (BHP). Let z G <9D, < r < r (D), and < p < 1. Ifu,v are 
nonnegative in lR d , regular L-harmonic in D(lB(z,r), vanish on D c (lB(z,r) 
and satisfy u(x ) = v(xq) for some x G D D B(z,pr) then 

C~ l v(x) <u(x) <Cv(x) , x G DnB{z,pr) , (63) 

wit/i C = C(d, a, 6,p, r (D)). 
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Proof. In view of Lemma [TS] we may assume that r is small. Let F = 
F(z,r/2) C B(z,r) be the C 1 ' 1 domain of Lemma [TJ localizing D at z. For 
x G F we have w(x) = / Pf(x, z)u(z) dz ~ where = / Pf{x, z)u(z) dz. 
Similarly v(x) ~ = J Pf(x, z)v(z) dz. Since £t(xo) = S(a:o), we have 
«(^o) ~ ^(^o)- By Lemma [31 u(x) ~ v(x), provided x G -D D B(z,r/8). We 
use Lemma [T51 for the full range x £ D C\ B(z,pr). □ 



5 Proof of Theorem [T] 

By fl52l) and (1291) we have the estimate 

G(x,y)<G(x,y) + d [ G j X ^ G{z ' y h b(z)\dz , x, y eD. (64) 

Jd A |y - z\ 

We consider 77 < 1, say 77 = 1/2. By Lemma [9] and the uniform integrability 
in Lemma HH (see (T49])) there is a constant r > so small that 

/ M z)\dz< n -, yeD, (65) 

J D r 5{z) A\y - z\ d 

and 

/ «%^I»WI*<1 veD. (66) 

,/ D r G(x,y)(5(z) A\y — z\) d 
Here _D r = {z G -D: £(2:) < r}. We denote 

p = [e Ar (D) Ar]/16, 

with e = e(d, a, b,2/n) of Lemma [131 see a l so Lemma [U 

To prove ([T]) we will consider x, y in a partition of D x D. We will also 
consider Q,R £ dD such that £(x) = |x — Q|, S(y) = \y — R\. 

I. First we suppose that S(y) > p/4. We denote 

• D 1 = {x G D: \y - x\ < p/8}, 

• D 2 = {xeD: 8{x) > p/8}, 

• D 3 = {x ED: 6{x) < p/8}. 
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a) Let x G D^. Denote B = B(y,p/4). By §M§ we have Ge(x,y) ~ \x — 
y\ a ~ d pa G(x,y). Lemma [T51 yields the lower bound in (1T1): 

y) > G B (x, y) « G B (x, y) « G(x, y) . 

By Lemma [7] we get the upper bound in (pQ). 

b) Let x G Z?2 \ -Di- Let xo G D be such that |xo — y\ = p/8. G(-,y) is 
L-harmonic in B(x,p/8) U B(xq, p/8). By Lemma [TBI a) and Lemma [2] we 
get G(x,y) « G(x ,y) « G(x Q ,y) « G(x,y). 

c) Let a; G -D3. Let xo, yo £ -D be collinear with x, Q and such that <5(xo) 
./',:;, — Q| = 5p/32 and <5(yo) = 1 3/0 — Q| — 7p/32 (consider an inner ball 
tangent at Q to see the situation). Let F = F(Q,p) be the approximating 
domain of Lemma [TJ The functions G(-, y ), Gp(-, Ho) an d G(-, y) are regular 
Inharmonic in B(Q,6p/32). By Lemma [T6| Lemma [TBI and Lemma [33 

G(x,y) ; G F (x,y ) G F (x,y ) ( ; G(x,y) 
G(x ,y) ~ G F (xo,yo) ~ Gf(zo,J/o) ~ G(x ,y) ' 

By a) and b), G(xo,y) ~ G(xo,y) and we obtain (JTJ in the considered case 
I, that is for <5(y) > p/4 and all x E D. 

Before we proceed to the next case we recall that G is non-symmetric. 

II. Suppose that 5(y) < p/4. The proof of ([TJ follows in 3 steps. 
Step 1. We will first prove that G(x,y) > cG(x,y), x G D. 
To this end we denote 

• F x = {x G D: \x-R\<p}, 

• F 2 = {x G D: 5(x) > p/4}, 

• F 3 = {x G D: S(x) < p/4}. 

d) Let x G F\. Consider F = F(R, 8p). Then S(x) = 5f(x) and 5(y) = 5p(y)- 
Consequently, G(x,y) ~ G F (x,y), see (1241) . As before we have 

G(x, y) > Gf(x, y) « Gf(x, y) « G(x, y) . 

e) Let x G F 2 \ i 7 !. Let 6 D be collinear with y, i? and such that 
5{xq) = \xq — R\ — p/2. We note that \x — y| > 3p/4 and \x — y| > p/4. By 
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Harnack inequalities and d), G(x,y) ~ G(x ,y) > cG(x ,y) ~ G(x,y). 

f) Let x G F 3 \ F\. Let 2 6 D be such that 5(zo) = \zq — R\ = p/3 and let 
xo G D be such that 5(xo) = |^o — Q| = p/4. We have |xo — j/| > p/2 and 
\zo — y\ > p/12. By Harnack inequalities and d), 

G(z , y) « G(*b, y) > cG(z , y) « G(x , y) . (67) 

G(-,2o) an d G(-,y) are regular L-harmonic in D fl B(Q,p/3) because |^o — 
Q| > hp/12 and |y — Q\ > p/2. By Lemma [T6| part I and Lemma El 

G(x,y) G(x,z ) G(x,z Q ) G(x,y) 



G(x ,y) G(x ,z ) G(x ,z ) G(x ,y) 
By this and ( l67l) . G(x,y) > cG(x,y). 

Step 2. We next prove the upper bound in (CD) for 5(x) > p/4. By part I, 

c^G{x, z) < G(x, z) < Cl G(x, z), zeD\D r . 

The constant C\ and other constants in what follows will only depend on d, 
a, (the suprema in the Kato condition for) b, r (D) and diam(D). 
By f|45|) and Lemma [TTj 

f G(x,z)G(z,y) 
A |y - z\ 

Therefore by (l64l). 

G(z,y)<G(z,y) + c lC / / I^MmW)]^ 

r G(x g )G(z >y) 
J DT 8{z) A \y — z\ 

<AG(x,y) + d[ G{x ' Z ] G } Z > V h b(z)\dz, (68) 
J D r \y - z\ A d(s) 

where A = 1 + c\dC\. By Lemma [7] and ( 1651) we obtain 

G(x,y)<AG(x,y) + C d [ ^ ~ f^ G{z ^ \b{z)\dz 

J D r 5{z) A \y - z\ 

< AG(x,y) + B(x) , (69) 
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where B(x) = r]Co5£> r {x) a d . We claim that for n = 0, 1, . . ., 

G(x, y)<A(l + r) + --- + r] n )G(x, y) + r) n B(x) . (70) 

This is proved by induction: we plug (JTD]) into fl68|) . and use fl65l) and (166|) . 
In consequence, 



G(x,j/)< 



.4 



1 — 77 



G(x,y) 



(71) 



Step 3. We will now prove the upper bound in ([I]) when 5(x) < p/4. We 
will consider Fi, F 3 and F from Step 1. If x G F 3 \ F\ than we use the same 
argument as in f), but this time all the terms in ( J67l) are comparable because 
of ( 17T1) . and we obtain (CQ). Finally, for x G i<i C F we have 



G?(a;,j/) = G F (x,y) 



Pf(x, z)G(z, y) dz . 



D\F 



By Lemma [TBI . G F (x,y) ~ G F (x,y). We already know that for z G -D \ iq, 
G(z,y) « G(z,y), and P F (x,z) « Pjrfo*) by (EUD. Thus, 

y) « G F (x, y) + / P F (x,z)G(z, y) dz = G(x,y) . 

J D\F 

The proof of Theorem [1] is complete. In passing we only note that (I38j) and 
([TJ, and [20] or [301 Theorem 22] yield sharp estimates of the Poisson kernel: 

P D (x,y) « P D (x,y) « 5 D (x) a/2 5 DC (y)^/ 2 [1 V My)P /2 |y - x\~ d . (72) 



6 Appendix 

Let function <fi be of class C ,x , i.e. satisfy 

|V0(£) - V0(y)| < 7}\x-y\, x,yeR d -\ 
for some 77 < 00. Let /(t) = 0((1 — t)y + £5), £ G [0, 1]. We have 



(z) -0(y) - V0(y) • 



(x - jjf) • (V0((1 - t)y + ix) - V0(y)) 



< 77|x — y| /2 . 



We will consider nonlinear transformations of lR d defined as follows, 
Tx = (x, x d + 4>{x)) , T~ x x = (x, x d - 4>(x)) . 



(73) 



(74) 



(75) 



27 



Proof of LemmaUl For x — (xi, . . . , x d _i,x d ) G R d we let x = (xi, . . . , x d -i), 
so that x = (x,x d ). The halfspace H = {x G M. d : x d > 0} is a C 1 ' 1 domain 
at each scale r > 0, and we can localize it at by 

K = {x eR d : < x d < r/2, \x\ < r/4 + v/(r/4) 2 - (x d -r/4) 2 } . 

Put differently, is defined by the conditions: \x\ < r/2 and 

r/4 - ^/(r/4) 2 - (|x| - r/4) 2 < x d < r/4 + ^(r/4) 2 - (|x| - r/4) 2 . 

We see that K is C 1 ' 1 at scale r/4. We consider the inner and outer balls of 
radius r/4 for K, tangent at Q G OK. Let 

B = {x G R d : 2(n, x - Q) > \x - Q\ 2 } 

be either one of them. Here n G M. d and |n| = r/4. 

Let D be C ' at a scale r > 0. It is well known that up to isometry D 
locally coincides with the image of the halfspace if by a transformation T of 
the form (f75l) . see [H Section 2]. Namely, by possibly changing coordinates, 
we may assume that G dD, satisfies ( T73|) . 0(0) = 0, V0(O) = 0, 

,~|2 

10(^)1 < < T ? provided |x| < r/2, 

r 4 

and 

{x E D : \x\ < r/2, \x d \ < r } = {x G R d : |x| < r/2, 0(5) < x d < r} . 

We also have t] < c/r in fl75|) . where c > 1 is an absolute constant. 

We define F = TK C -D. We see that F locally coincides with D, and 
F C {x : \x\ < r/2, —r/4 < x d < 3r/4}, hence diam(F) < 2r. We claim 
that F is C 1 ' 1 . The claim will follow from considering the image of B by T. 
Let f = V0(Q). We note that |f | < < r/r/2 < c/2. We define 

x ^ Sx = T(Q + x) - TQ = (x, x d + [0(Q + x) - 0(Q)]) , 

5' _1 x = (x,x d - [cj)(Q + x)-(j)(Q)}), Lx = (x,x d + £-x), L~ l x = (x,x d -^-x). 
We will use L^ 1 as a linear approximation of S*" 1 at x — 0. Let be the 
transpose of L -1 . We have \Lx\ < \x\ + |^||x| < |x|(2 + c)/2. The same is 
true of L _1 and L~ l * . We note that 

TB - TQ = S(B - Q) = {x G R d : 2(n, S ,_1 a:) > l^af } , 
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and x = is on the boundary of the set. Consider the ball B' = {x : 
2k (L _1 *n, x) > \x\ 2 }. Note that TQ G d{B' + TQ). We will verify our claim 
on F by proving that if 2k = c _1 A (1 + 3c/4)" 2 , then B' C 5(5 - Q), or 
i?' + TQ C T£>. To this end we note that 

2K(L~ u n,x) > \x\ 2 (76) 

implies that \x\ < 2/t|L~ 1 *n| < r/2. For such (small) x, by (|74p . we obtain 

l^-^l < |x| + + r/|x| 2 /2 < |x|(l + 3c/4) . (77) 

Similarly, 

|5 _1 a; - L~ l x\ = \<p(Q + x) - <f>[x) -£-x\< r]\x\ 2 /2 . 
Now, (J75]) yields 2(n,L~ 1 x) > |x| 2 //t, hence 

2(n, S^x) > 2(n, S~ l x) - 2(n, L -1 x) + -Id 2 > -b| 2 - 2|nb|x| 2 /2 

K K 

> —\x\ 2 > (2k)- 1 (1 + 3c/4y 2 \S~ 1 x\ 2 > IS^xl 2 . 
2k 

The proof is complete. □ 
Proof of LemmalE By [S] for s > 0, x G M d and G C c °°((0, oo) x R d ), 

p(u — s, x, z) {d u + A"/ 2 + b(z) ■ V z ) 4>(u, z)dz du = —<p(s, x) . (78) 

We note that the above integral is absolutely convergent. Indeed, 

| (d u + A"/ 2 + b(z) • V 2 ) cj>(u, z)\ < c(l + |6(*)|) , 
and if z) = for u > M, then by (133]) and the remark following (1301) . 

np(n - s, x, z) | (d u + A q/2 + b(z) ■ V z ) 0(n, z) | dz du 
— c i / / p(u — s,x, z)(l + \b(z)\) dzdu 

Js JR d 

<c 2 [ (\x - z\ a ~ d A\x- z\- a ~ d )(l + |6(z)|) dz < oo . 
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Furthermore, for x G s G R by and (I78p we obtain 
y y P£>(w - s > ^) (<9u + A" /2 + 6(2;) • V 2 ) (j)(u, z) dz du 

POO P 

= / p(u-s,x,z) (d u + A* /2 + b(z) ■ V 2 )0(m, z)dzdu 

Js JR d 

roc p 

-W p{u-s- r D ,X TD ,z) (8 U + Af + b(z) ■ Vj </>(u, z) dzdu 

Js+t d JR d 

= -<f>(8,x)+W<j>(8 + T D ,X TD ) (79) 

□ 

Proof of Lemma\^ Let s = 0, x G -D, G C£°((0, 00) x M d ) and assume that 
supp G (0, 00) x (U) c . By ([7SD and © we obtain 



E x 0(r D ,X TD )= / / / p D (u,x,y)(j)(u,z)u(z - y)dzdydu . 

JO J D JD C 



□ 



Proof of Lemma\B Let <p G C~(D) and let x(u) G C C °°(IR) be such that 
x(u) = 1 for u G (—1, 1). For rt = 1, 2, ... we define (f) n (u, z) = ip(x)x(u/n), 
and we have 

I (d u + A a/2 + &(*) • V,) n («, z)| < c(l + |6(z)|) 
with c independent of n. By ( 1361) . Lemma [7] and dominated convergence, 
- tp(x) = lim -0 n (O,x) 

n— >oo 
poo r 

lim / / po(u, x, z) (d u + A"/ 2 + • V z ) 4> n (u, z) dzdu 
n ^°°Jo Jd 

r"0O /• 

p D (u, x, z) (A a/2 (p(z) + b(z) ■ V(f(z)) dz du 
G D (x, z) (A a / 2 cp(z) + b{z) ■ Vcp(z)) dz . 

□ 




L) 
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